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1 Introduction 

The study of cold quark matter has been an interesting topic of research in recent 
years. In an astrophysical context, quark matter may be used to explain the mea- 
sureable properties of pulsar-like compact stars. 

It has recently been proposed that realistic quark matter in compact stars could 
exist in a solid state O [3], either as a super-solid or as a normal solid j7]. The 
basic conjecture of normal solid quark matter is that de-confined quarks tend to 
form quark-clusters when the temperature and density are relatively low. Below 
some critical temperature, these clusters could be in periodic lattices immersed in a 
degenerate, extremely relativistic, electron gas. Note that even though quark matter 
is usually described as weakly coupled, the interaction between quarks and gluons 
in a quark-gluon plasma is still very strong [5]. It is this strong coupling that could 
cause the quarks to cluster and form a solid-like material. 

In this paper, we argue that comtemporary observations of thermal X-ray emitting 
pulsars are consistent with the assumption that these sources are in fact Solid Quark 
Stars (SQSs). In the successive Sect., the thermal X-ray observations on isolated 
pulsars are exhibited according to the collation in this work. In Sect. 3, the SQS pul- 
sar model is explored to interpret these observations. Conclusions and corresponding 
discussions are presented in Sect. 4 



2 The observations 

We collate the thermal X-ray observations on 29 isolated pulsars (see Table 1 in [9]). 
The black body parameters summarized comprise of pulsar surface temperature com- 
ponents T s °^yJ2l, emission size components -R^ 2 ' a § es anc ^ bolometric X-ray luminosities 

1 "oo" — values measured at distant. 
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Figure 1: Functions L™ ol (E) (top panel) and t(E) (bottom panel) of active pulsar 
candidates. In the top panel, the fits are carried out both for Group A (including top 
17 pulsars listed in Table 1 in [9]; they are marked by dark points and their numbers) 
and Group B (including the members in Group A and the other 10 sources whose 
upper limits on the bolometric luminosities could be defined observationally; the 10 
sources are marked by grey points and letters). Red lines are the fitted results for 
Group A, while blue lines are those for Group B. In both groups, solid lines provide 
the best fits to the data, while the dashed lines give the fits by freezing pi at 1. We 
note that the 10 sources with upper limits on their bolometric luminosities defined are 
taken from p, and they are a. B1509-58, b. B1951+32, c. B1046-58, d. B1259-63, 
e. B1800-21, f. B1929+10, g. B0540-69, h. B0950+08, i. B0355+54, j. B0823+26. 
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0.4561 ±0.2315 
(1) 

u.t>9±8_ 2046 
(1) 


16.20 ± 8.30 

n 9 on+0.495 

O.^OU_q 

11 66+ 73 ' 100 
11 '™-7 .2990 

-2.896 ± 0.403 


0.7487 
0.7730 


0.3277(14) 
0.8607(15) 
0.6081(24) 
0.9964(25) 



a Values in the parenthesis are frozen during the fits. The errors are in 95% confidence level. 
b Correlation coefficient of the data between bolometric luminosity and E. 
c Reduced x 2 i or \ 2 per degree of freedom (d.o.f). 

Table 1: Fitting parameters for L™ ol — E data, which has been shown in Fig. [1] (top 
panel) . 



-^boi- Spin would be a significant source of energy for pulsars. Potential relation be- 
tween pulsar X-ray bolometric luminosity L^ ol and its spin energy loss rate E could 
exist. Such a study was done and the results are exhibited in Fig. [I] and Table [TJ 
Using the function of 

\ogLZ l =pl-\ogE + p2, (1) 

L^ ol — E data were fitted and the results show that the best fit demonstrates a possible 
1/2-law between these two quantities, or 

LZi(E) = C&l\ (2) 

where the coefficient C = 10 p2 is in the unit of erg 1 / 2 s _1 . However, when the param- 
eter pi is set to 1, a linear-law of 

LZm = r ] -E, (3) 
is presented, and this coefficient rj is 10~ 3 approximately. 

3 The SQS pulsar model 

3.1 X-ray cooling and accreting pulsars 

The state of a pulsar would be manifested in multiple wave bands and by different 
photon components. Top 17 pulsars have significant observable X-ray non-thermal 
components, which could originate from energetic magnetospheres. Such pulsars are 
considered to be undergoing cooling processes as proposing by SQSs. The rest sources 
are dominated by thermal X-ray components observationally, and show little magne- 
tospheric activities. These pulsars may be accreting to the surrounding interstellar 
medium to generate visible X-ray luminosities, if they are in fact SQSs. 
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3.2 Cooling of SQSs 
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Figure 2: Left panel: Cooling curves for SQSs, if the 1/2-law between the bolometric 
luminosity and the spin energy loss rate holds. Right panel: Corresponding tempera- 
ture differences between the hot and warm components of SQSs, or AT = T£° — T s °°. 
The parameters in both panels: M = 0.1M Q , C = 10 16 (solid lines); M = 1.0M Q , 
C = 10 16 (dashed lines); M = l.OM , C = 10 15 (dotted lines); M = O.O1M , 
C = 10 15 (dash-dot lines) (C is in the unit of erg 1 / 2 s™ 1 / 2 ). For two curves with 
same M and C, the upper one corresponds to an initial spin of 10 ms, while that of 
the lower one is 100 ms. Noting that the errors on the surface temperatures of PSRs 
J0205+6449 (No. 3) and J2043+2740 (No. 17) are not provided by the authors of the 
references, we then conservatively adopt them as deviating from the central values by 
a factor of 2. 



Cooling of a quark star, from the initial birth, can be quite complicated to model. 
Two facts induce this result. Firstly, recent experiments on collisions between heavy 
relativistic ions (e.g. RHIC) manifested that the potential quark-gluon plasma in the 
fire ball would be strongly coupled even in the temperatures of a few factors of the 
critical temperature ~170 MeV [5]. However, and secondly, the initial temperature 
at the birth of a pulsar could only be a few factors of ~ 10 MeV. Hence, if the 
natal pulsar is a quark star, the stellar interior quark-gluon plasma would be coupled 
strongly, so that current physics could not model either in an analytic or a numerical 
way. 

Nevertheless, a phenomenological study on this problem could be carried out, and 
Fig. 1 in [9j illustrates proposed cooling processes for a quark star. It could be quite 
uncertain that at which temperatures the phase transitions would occur, and how 
long the star would stay in Stages 1 and 2. However, in consequence of the strong 
coupling in quark-gluon plasma, Stages 1 and 2 would not be long, and the observed 
pulsars (if they are quark stars) would be in the solid phase (Stage 3), because of the 
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Figure 3: Left panel: Cooling curves for SQSs, if the linear-law holds. Right panel: 
Temperature differences for this case. The parameters: M = l.OM , r] = 0.01 (solid 
lines); M = l.OM , rj = 0.001 (dashed lines); M = O.1M , r] = 0.1 (dash-dot lines); 
M = 0.01M Q , T] = 0.1 (dotted lines). As in Fig. [21 for two curves with same M and 
r), the upper one corresponds to an initial spin of 10 ms, while the lower one 100 ms. 



measured low temperatures. The calculations below concentrate on the late Stage 3, 
where the surface temperature is in the range of the observed 10 7 — 10 5 K. The loss of 
the thermal energy of a SQS in the late 3rd stage could dominantly induced by the 
emission of photons rather than neutrinos, as a result of the low neutrino emissivities 
in such low temperatures. 

In general, the cooling process would be described by the equation of 

— C v — — + L SH = L ho \, (4) 
at 

with assuming that the stellar volume is invariant and the internal energy of the star 
is only the function of the stellar temperature T s . C v is the stellar heat capacity. Lsh 
is the luminosity of the probable heating process, and Lboi is the photon luminosity 
of the star. The stellar heat capacity was found to be tiny, so that the internal 
energy would be far from sustaining long-term X-ray emission. We note here that 
for a 1.4M -SQS, it would cool down to 1 K from 10 11 K in ~0.2 s. The detailed 
exhibition of the calculation of this stellar residual internal energy as well as heat 
capacity can be found in Appendix A. 

Heating mechanism would take important effect during the cooling of SQSs. In 
such a way, the thermal evolution of a SQS could not be a self-sustaining cooling 
process, but an external-supplying heating process. The bombardment by backflowing 
particles at the poles could be the prime heating mechanism for SQSs, which would 
intrinsically be spin-powered. Combining the negligible residual internal energy of 
a SQS and the observed L^ ol — E relations. The heating luminosity Lsh could also 
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follow that 1/2-law or linear-law. The equation would hence yields 



CE 1 ' 2 



4nR 2 aT s A + ^r\aT^ 



(5) 



or 



rjE 



47rR 2 aT* + 4iir 2 aT A 

b P P 



(6) 



The two terms on the right side of the equal sign are the stellar primary black body 
emission component and that for a probable hot pole component respectively. The 
relation between T p and T s could be found by 



where H is the heat current flowing from the poles to the bulk of the star. The stellar 
heat conductivity n c mainly contributed by free electrons refers to the corresponding 
results for solid metal [2]. Figs. [2] and [3] exhibit the computed cooling curves for the 
above two cases respectively. In these calculations, SQSs were assumed to slow down 
as orthogonal rotators that lose their spin energy by magnetic dipole radiation. The 
average densities of the SQSs were adopted to be 3 times nuclear saturation density. 
Note that the results will not vary significantly when the density values are in the 
range of 3-5 nuclear saturation density, which is thought to be representative values 
for the solid quark matter. The spin-powered thermal emission of SQSs enable an 
estimate on pulsar moments of inertia [9]. 



Dim X-ray isolated neutron stars (XDINs) and compact central objects (CCOs) 
(Sources No. 18-29 and 4) are thought to be inactive pulsar candidates with inert 
magnetospheric manifestations. These sources may have excessive X-ray luminosities 
to their spin energy loss rates. Their energy origin has been keeping disputable, but 
accretion origin would be a candidate. The relative high magnetic fields of these 
sources may drive their surrounding shells of matter (in large radii, usually larger 
than the corotation radii) to rotate together with them. Under this propeller regime, 
nearly all of the matter would be expelled outward, but a small fraction may diffuse 
starward and fall onto the stellar surfaces eventually. The release of the gravitational 
energy as well as the latent heat inducing by the burning from the baryonic mat- 
ter phase to the three-flavor quark matter phase could sustain long-term soft X-ray 
radiations of such pulsars. In this scenario, the stellar luminosity would then be 
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3.3 Accreting SQSs 
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For the calculated accreting parameters for SQSs in this propeller regime and the 
resulting luminosities, please see [9]. Note here that for No. 4 RX J0822-4300, a 
CCO, we treated it both as a cooling SQS and a accreting SQS, as the result of the 
disputable structure around it pUl H] . 

4 Conclusions 

We collate the thermal observations of 29 X-ray isolated pulsars and, in the SQS 
regime, for the magnetospherically active pulsar candidates, establish their cooling 
processes (Figs. |2]and[3]), while for the magnetospherically inactive or dead pulsar 
candidates, interpret the X-ray luminosities under the accretion scenario. SQSs, 
because of the possibility of being low-mass, could provide an approach to understand 
the small black body emission sizes. We note that for SQSs with mass of O.O1M , 
O.1M and 1M , their radii are ~1.8, ~3.8 and ~8.3 km respectively. On the other 
hand, a linkage between pulsar rotational kinetic energy loss rates and bolometric 
X-ray luminosities is explored by SQSs. 

We hence conclude that the phenomenological SQS pulsar model could not be 
ruled out by the thermal observations on X-ray isolated pulsars, though a full depic- 
tion on the formation and the thermal evolution for a quark star in all stages could 
hardly be given nowadays, as the lack of the physics in some extreme conditions (as 
has been discussed in Sect. 3.2). SQSs have significant distinguishable interiors with 
neutron stars, but the structures of the magnetosphere between these two pulsar mod- 
els could be similar. Therefore, SQSs and neutron stars would have similar heating 
mechanisms, such as the bombardment by backflowing particles and the accretion to 
surrounding medium. However, a full comparison between SQSs and neutron stars 
including their cooling processes as well as heating mechanisms has beyond the scope 
of this paper, and should be interesting and necessary in the future study. 

A Residual internal energy of a SQS — fast cool- 
ing? 

Provided that the volume of a SQS is a constant, then the stellar residual thermal 
energy £7sqs would only be the function of the stellar temperature T s , 



where T s is the value in the star's local reference frame, and C v is the heat capacity of 
the star. Using Debye elastic medium theory, the characteristic of the heat capacity 




(9) 
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of solid quark matter would be evaluated by Debye temperature 6d, or 

e D = hu D /k B , (io) 

where h is reduced Planck constant, and &b is Boltzmann constant, ujd is Debye 
cut-off frequency (i.e. the maximum frequency of the wave that could propagate in a 
medium), which equals Debye wave number kr> = (67r 2 n c ) 1,/3 times the average sound 
speed in the medium, i.e. ujd = ^dc s , where n c is the number density of classical 
particles, or quark clusters for solid quark matter. For a SQS, the average sound 
speed could be the light speed approximately. A linear equation of state, extended 
to be used for a quark star, indicates that the pressure p ~ pc 2 , where p is the 
mass density of a quark star and c is the speed of light. So an estimate would be 
c s = a/ dp/ dp ~ ypjp ~ a/ pc 2 / p = c. n c = 3en /A, where e denotes the baryon 
number density of solid quark matter in the unit of no- no = 0.17 fm -3 is the nuclear 
saturation baryon number density. We consider e = 3 — 5 could be the representative 
values for a SQS. A is the number of valence quarks in a quark cluster. We may expect 
that A could be in the order of 10, since A = 18 if quark-a-like clusters are formed [6], 
and A could even be conjectured to be in the order of > 10 2 . Debye temperature 6d 
of a SQS could then be of order ~ 10 12 K, even much higher than the temperature 
when a quark star is born. Hence, the heat capacity in the low-temperature limit or 
the temperature-cube law, would be applicable for the SQS stage (Stage 3), i.e. 

c v = ^b(^) 3 , (11) 
5 D 

where c v is the heat capacity per classical particle (or quark cluster when refers to 
the solid quark matter). Thus C v = N ■ c v , where N is the total number of clusters 
in a star. 

By referring to equation (jUJ), the residual internal energy for a 1.4M Q -SQS would 
thus be Usqs ~ 10 4 T S 4 erg, while the corresponding value for a typical neutron star 
in such mass is U^s ~ 10 30 T 2 erg [S]. Therefore, the residual internal energy of a 
SQS will always be much less than that of a normal neutron star; even when the 
temperature is as high as ~ 10 11 K, the residual thermal energy of a SQS is only one 
ten thousandth of that of a neutron star. So if the cooling process of a 1.4M -SQS is 
only powered by its residual internal energy, i.e. 

- C v ^ = 4nR 2 aT* + ^i? 3 10 24 T s 6 9 (12) 

(R is the local stellar radius, a is Stefan-Boltzmann constant), then its cooling 
timescale from 10 11 K to an extremely low final temperature 1 K is only ~ 0.2 s. 
Noting that as a result of the concerning of the high initial temperature here, neu- 
trino luminosity for the SQS is also considered as the second term on the right side 
of the equal sign in equation ffT2]) . 
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